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We apply a model for sand dunes to calculate formation of dunes on Mars under the present 
Martian atmospheric conditions. We find that different dune shapes as those imaged by Mars 
Global Surveyor could have been formed by the action of sand-moving winds occuring on today's 
Mars. Our calculations show, however, that Martian dunes could be only formed due to the higher 
efficiency of Martian winds in carrying grains into saltation. The model equations are solved to 
study saltation transport under different atmospheric conditions valid for Mars. We obtain an 
estimate for the wind speed and migration velocity of barchan dunes at different places on Mars. 
From comparison with the shape of bimodal sand dunes, we find an estimate for the timescale of 
the changes in Martian wind regimes. 
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I. INTRODUCTION 

Sand dunes are ubiquitous on Mars and provide evi- 
dence that sand-moving winds have once occured on the 
red planet. However, Martian dunes do not appear to be 
moving these days, and thus many authors suggested that 
dunes on Mars have been formed in the past, when the 
climate of Mars was much more earthlike [1| . On a planet 
where the present atmospheric density is almost hundred 
times lower than the Earth's, only winds one order of 
magnitude stronger than on Earth are able to transport 
sand and form dunes [3, . Such strong winds occur in 
fact only occasionaly on Mars, during the strongest dust 
storms [3, Q . Only once during the Viking mission, which 
operated for about six terrestrial years, could very little 
changes on Martian soils be detected after such a storm 
0. Moreover, calculations showed that the sand-moving 
winds that occurred at the Viking landing site lasted for 
not more than a few tens of seconds 0, 0] . Could Martian 
dunes have been formed by the action of such rarely oc- 
curing winds under the atmosphere of the present Mars? 
This is the question that motivates the present work. 

The sand that forms dunes is transported by the wind 
through saltation, which consists of grains travelling in a 
sequence of ballistic trajectories and producing a splash 
of new ejected grains when impacting onto the ground 
Q. Martian saltation has been studied in wind tunnel 
experiments and also in numerical simulations 0, Q . Due 
to the thinner atmosphere of Mars and owing to the lower 
gravity g = 3.71 m/s^, which is nearly 1/3 of the Earth's 
gravity, saltating grains on Mars travel longer and higher 
than their terrestrial counterparts Q. Moreover, Martian 
grains also saltate with faster velocities than grains on 
our planet. As consequence, the grain-bed collisions or 
splash events on Mars are expected to be much larger 
than on Earth, due to the larger momentum transfered 
by the impacting grains to the sand bed [10]. 

What is not known is whether such highly energetic 
saltation events have been responsible for the formation 



of the enormous dunes observed in the images of Mars. 
In order to understand this, it is necessary to investigate 
sand transport at length scales comparable to the scale of 
dunes. Once saltation starts, the wind transfers momen- 
tum to accelerate the grains. Thus, the momentum of the 
air decreases as the flux of saltating particles increases 
("feedback effect" [HI). After a distance which is called 
"saturation length" , the wind strength is just sufficient 
to sustain saltation, and the sand flux achieves satura- 
tion. In this manner, dunes that have length smaller than 
the saturation length will be continuously eroded due to 
increase of the sand flux and will disappear. In other 
words, the existence of a minimal dune size is related 
to the phenomenon of flux saturation, which could not 
be investigated from wind tunnel simulations of Martian 
saltation [l^]. While the first wind tunnel simulat- 
ing Martian conditions is a few meters long, the smallest 
dunes on Mars have length of the order of hundred meters 

(fig. m- 

Recently, a successful modelling of the formation of 
sand dunes, which encompasses the main processes of 
saltation and accounts for flux saturation and the exis- 
tence of a minimal dune size, has been achieved [H, [11] . 
This model consists of a system of continuum equations 
in two space dimensions which reproduce the shape of 
terrestrial dunes, the wind profile and the sand flux and 
provide excellent quantitative agreement with measure- 
ments fl^. The dune model, which has been applied to 
study the interaction of dunes in a field [IB] and the for- 
mation of parabolic dunes in the presence of vegetation 
[l^, has become a powerful tool in the investigation of 
the large timescale processes involved in the formation of 
desert and coastal landscapes. 

In the present work, we apply the dune model to in- 
vestigate whether dunes could be formed on the present 
Mars. Our aim to reproduce the shape of Martian dunes 
using the present Martian atmospheric conditions. The 
dune model has parameters of atmosphere, wind and 
sand, many of which are known for Mars and can be 
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FIG. 1: Mars Global Surveyor (MGS), Mars Orbiter Camera (MOC) images of sand dunes on Mars (courtesy of 
NASA/JPL/MSSS). From the left to the right: a. Barchan dunes at Arkhangelsky crater, near 41.0°S, 25.0° W; b. north polar 
dunes near 77.6°N, 103.6° W; bimodal sand dunes near c. 48.6°S, 25.5° W; d. 49.6°S, 352.9° W, and e. 76.4°N, 272.9°W. 



therefore used to solve the model equations. While most 
of the quantities controlling saltation can be calculated 
from the atmospheric density pfluid, gravity g, air viscos- 
ity rj and from the grain diameter d, there is one unknown 
quantity, which is related to the intensity of the Martian 
splash and must be determined from simulations. More- 
over, the wind velocity that formed Martian dunes is also 
a parameter: it must be estimated from comparison with 
the shape of the dunes. 

This paper is organized as follows. In the next Sec- 
tion we describe the dune model. The main equations 
are presented, as well as the relations used to calculate 
the microscopic quantities governing saltation on Mars. 
In Section IIIII we calculate the average grain velocities, 
mean saltation height and sand flux on Mars from the 
known parameters, and the results are compared with 
wind tunnel predictions and with calculation results for 
saltation transport on Earth. Dunes on Mars are calcu- 
lated in Section IIVI We begin with the simplest dune 
form, which is the barchan dune. First, we study the 
shape of the barchan dunes in the Arkhangelsky crater. 
We find an equation for the rate at which grains enter 
saltation, which can be used in the calculations of dunes 
under different atmospheric conditions. We then esti- 
mate the wind velocity on Mars and predict the migra- 
tion velocity of Martian barchans. Next, we study the 
shape of Martian bimodal sand dunes and find an es- 
timate for the timescale of changes in wind regimes on 
Mars. Conclusions are presented in Section [V] 



II. THE DUNE MODEL 

The dune model combines an analytical description 
of the average turbulent wind velocity field above the 
dune with a continuum saltation model which allows for 
saturation transients in the sand fiux. In the model. 



avalanches on the slip face and fiow separation at the 
dune lee are also taken into account. Here we give a brief 
presentation of the model and refer to Sauermann et al. 
[ist and Schwammle and Herrmann for the extensive 
derivation of the saltation transport equations. 



A. Wind 

In the turbulent boundary layer, where sand transport 
takes place, the velocity of the wind u{z) increases loga- 
rithmically with height z above the fiat ground: 



u{z) = — m — , 
K zo 



(1) 



where n = 0.4 is the von Karman constant, m, is the wind 
shear velocity, which is used to define the shear stress 
T = PfluidU*, and Zo is the aerodynamic roughness, 
and zq are two independent variables that can be deter- 
mined from measurements of the wind velocity at differ- 
ent heights, as done for instance in the Mars Pathfinder 
Lander Wind Sock Experiment [18,]. The aerodynamic 
roughness zq is larger than the surface roughness of the 
undisturbed sand bed, zg™'^, which is of the order of a 
few tens of microns and is due to the microscopic fiuctu- 
ations of the sand bed when the grains are at rest . A 
value of zq close to 1.0 mm has been often reported from 
measurements of terrestrial saltation on a sand bed , 
while on Mars zq is larger, around 1.0 cm jT8| . 

A dune or a smooth hill introduces a perturbation in 
the shear stress whose Fourier-transformed components 
are calculated using the algorithm of Weng et al. : 



2 h(^kx^ ^y)^x 
\k\UHl) 



2ln{C\kx\)+4:e + l+i sign(fc^ 



ln(Z/zo) 



(2) 
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and 



2 h(^kxj ky)kxky 
\k\U^{l) ' 



(3) 



where the coordinate axes x and y are parallel, respec- 
tively, perpendicular to the wind direction, kx and ky 

are wave numbers, \k\ — ^k^. + fc^ and e — 0.577216 

(Euler's constant). £ is the horizontal distance between 
the position of maximum height, -ffmax, and the position 
of the windward side where the height is i?niax/2 [20j . 
U(V) — u{l)/u{hm) is the undisturbed wind velocity at 
height l — 2K^C/h\l/zQ normalized by the velocity at the 
reference height h-a^ = LI a/Io r Cj jzn, which separates the 
middle and upper flow layers [20| ■ The shear stress in the 
direction i {i = x, y) is then given hy Ti — i [to(1 + f^)], 
where tq is the undisturbed shear stress over the flat 
ground. 

Sand transport occurs if exceeds a threshold veloc- 
ity for entrainment, u*ft, which depends on d, g, pfiuid, 
on the grain density Pgrain and also on the packing of 
the grains [2l|. Indeed, the wind velocity may even de- 
crease to values lower than u*ft, and still saltation can 
be sustained, once initiated. This is because the splash 
is the most important mechanism of sand entrainment 
during saltation Q . The wind strength, however, cannot 
be lower than the impact threshold velocity u*t, which 
defines the threshold shear stress rt — PfiuidU*t and is 
around 80% of M*ft. Saltation ceases if < u^^t, and 
therefore the impact threshold velocity is the essential 
threshold velocity for aeolian sand transport. 



B. Continuum saltation model 

The wind shear velocity computed above is used to 
calculate the sand flux with the model derived in Sauer- 
mann et al. [l3| . The fundamental idea of the model is 
to consider the bed-load as a thin fluid-like granular layer 
on top of an immobile sand bed. 

The sand bed represents an open system which can ex- 
change grains with the moving saltation layer, for which 
the erosion rate r(a;, y) at any position {x, y) represents 
a source term. Mass conservation yields that the local 
change in the flux balances the erosion rate: 



v-[pv) = r(x,?/). 



(4) 



where p(a;, y) is the density of grains in the saltation layer, 
and v{x,y) is the average local velocity of the saltating 
grains, whereas stationary condition has been assumed 
{d/dt = 0) since the time scale of the surface evolution is 
several orders of magnitude larger than the typical values 
of transient time of the saltation flux. The erosion rate is 
the difference between the vertical flux of ejected grains 
and the vertical flux cj) of grains impacting onto the bed: 



(5) 



where n is the average number of splashed grains. The 
vertical flux (p is defined as (p — p\v\/£, where £ is the av- 
erage saltation length. Due to the multiplicative process 
inherent in the splash events, the number of saltating 
grains first increases exponentially. However, the "feed- 
back effect" leads to a decrease of the momentum of the 
air as the population of grains increases. At saturation, 
the air shear stress Ta decreases to the threshold rt , while 
the flux of grains increases to its maximum value. The 
number of ejecta compensates the number of impacting 
grains: n = 1 at saturation. In this manner, we write n 
as a function n(ra/rt) with n(l) = 1. Expansion of n into 
a Taylor series up to the first order term at the threshold 
yields 



(6) 



n= l + 7( -^-1 ), 



where 7 — dn/d(Ta/Tt), the entrainment rate of grains 
into saltation, determines how fast the system reaches 
saturation J^] . Inserting eq. ^ into eq. ([5|) , and substi- 
tuting this into eq. ([4]), we obtain a differential equation 
for the sand flux. 

However, the air shear stress Ta within the saltation 
layer is written as Ta = t — Tg, where Tg is the con- 
tribution of the grains to the total shear stress at the 
ground [ll]. The "grain born" shear stress is calculated 
as Tg = 0Awhor — p\v\Avhor/£, whcrc Auho 



^hor 



"hor 



gives the difference between the horizontal velocities (in 
the direction of the flow) of the grains at the moment of 
impact, f™^, and at the moment of ejection, w^o" . The 
equation of mass conservation can be written, then, in 
the following manner: 



, ^, p\v\ _ T — 
£ Tt 



1 - p\v\ 



n 



(7) 



On the other hand, the mean saltation length £ is defined 
in terms of the grain velocity \v\ and of the initial vertical 
velocity, vf'^ [l3|: £ = vf'^{2\v\/g), where is related 
to the gain in horizontal velocity of the grains, Awhor, 
through an effective restitution coefficient for the grain- 
bed interaction, a Il3ll. which is defined as 



,,cje 



Awho 



(8) 



In this manner, the mean saltation length is written as 
.2|iT| a 2\v\ 1 



9 



Avhoi g 



2\v\^a 



9 



where 



Awh, 



(9) 



(10) 



is the constant of proportionality between the average 
grain velocity \v\ and the difference between impact and 
ejection velocity of the grains, Awhor- 
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The velocity of the saltating grains, v, is determined 
from the balance between three forces: (i) the drag force 
acting on the grains; (ii) the bed friction which yields 
the loss of momentum when the grains impact onto the 
ground, and (iii) the downhill force, which acts on the 
saltation layer in the presence of bed slopes. To calcu- 
late the grain velocity, we need to take into account the 
modification of the air flow due to the presence of the 
saltating grains. However, the model equations do not 
account for the complex velocity distribution within the 
saltation layer. Instead, a reference height zi is taken, 
between the ground at the roughness height Zq™*^ and 
the mean saltation height Zm, at^which the "effective" 
wind velocity, UeS, is calculated |l3]. 

Finally, a useful approximation is employed, which 
simplifies the equations in significant manner and leads 
to only a negligible error [l^. In the model, the velocity 
UcS is calculated in the steady state, i.e. it is the reduced 
wind velocity within the saltation layer at saturation. In 
geomorphological applications, the sand flux is nearly ev- 
erywhere saturated, with exception of those places where 
external variables change discontinuously, as for instance 
at a flow-separation, which occurs at the dune brink, or 
at a phase boundary bedrock/sand which occurs at the 
windward foot of a barchan dune. Therefore, we can re- 
place the density p which appears in the expression for 
the grain born shear stress Tg with the saturated density 
Ps = {t — Tt)£/Az;hor .133- The following expression is 
obtained, in this manner, for UcS- 



state. 



UeS 



— < In 

1 .vSand 




-1-1 




The grain velocity, v, is, next, calculated numerically 
from the equation [l3|: 



3 Pfluid Cd 

4 Pgrain d 



(Uoft - v)\Ue« - V\ 



2a\v\ 



gVh = 0, (12) 



where C'd is the drag coefficient, and Ues is calculated 
with eq. (Ilip . In this manner, the grain velocity obtained 
from eq. (fT2l) is in fact the average grain velocity at the 
steady state, Vg, since Ucs is the reduced wind velocity 
after flux saturation has been achieved fl3|. Moreover, 
the mean saltation length (eq. ^) is also computed 
using \v\ = \vs\ in the stationary condition where the 
sand flux is in equilibrium. The average grain velocity 
Vs can be now substituted into eq. ([7]), which is then 
written in terms of the sand flux 



Q = PVs, 

and the saturated sand flux qs — Psl^^l, or 



'2a\vs\ 
Qs = [r-n) 



2a\v. 
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1 



(13) 



(14) 



The resulting equation for the sand flux is a differential 
equation that contains the saturated flux qs at the steady 



V-q 



(15) 



where £s — [^/7]''"t/(''" — Tt) is the saturation length, which 
contains the information of the saturation transient of the 
sand flux. Using eq. is may be written as 
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(u^/u^t) -1 
where we deflned 

"f — rj 



2|^4| a/g 
(u*/u*t)^ 



1 



At;ho 



dn 



(16) 



(17) 



C. Surface evolution 

The change in the surface is computed using the flux 
calculated with eq. (fT5|) . The surface is eroded wherever 
the sand flux increases in the direction of wind flow, and 
sand deposition takes place if the flux decreases. The 
time evolution of the topography h{x,y,t) is given by 
the mass conservation equation: 



dh 
'dt 



1 



Ps 



-V-g, 



(18) 



where psand = 0.62pgi.ain is the mean density of the im- 
mobile sand which constitutes the sand bed [l3j- If 
sand deposition leads to slopes that locally exceed the 
angle of repose 9,- ~ 34°, the unstable surface relaxes 
through avalanches in the direction of the steepest de- 
scent. Avalanches are assumed to be instantaneous since 
their time scale is negligible in comparison with the time 
scale of the dune motion. At the brink of the dune, which 
represents a sharp edge, there occurs flow separation. In 
the model, the separation streamlines are introduced at 
the dune lee as described in details in Kroy et al. p^ . 
Each streamline is fltted by a third order polynomial con- 
necting the brink with the ground at the reattachment 
point ^12,] , and deflning the "separation bubble" , in which 
the wind and the flux are set to zero. 

The dune model can be sketched as follows: 

1. the shear stress over the surface is calculated with 
the algorithm of Weng et al. using eqs. ^ 
and (131); 

2. from the shear stress, the sand flux is calculated 
using eq. (|15p . where the saturation length £s and 
the saturated sand flux qs are calculated from ex- 
pressions (HH) and (HU, respectively; 



3. the change in the surface height is computed from 
mass conservation (eq. (IT51) ') using the calculated 
sand flux; and 
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4. avalanches occur wherever the incHnation exceeds 
34°, then the shp face is formed and the separation 
streamhnes are introduced as described in Kroy et 
al. [13. 

Calculations are performed using open boundaries with 
a constant influx of sand, gin, at the inlet. The influx is 
interpreted as the average interdune flux in a dune fleld, 
which is typically between 10 and 40% of the maximum 
flux qs (22.] . and is considered, for simplicity, homoge- 
neous along the y axis (perpendicular to sand transport). 
The model is evaluated by performing steps 1) through 
4) computationally in a cyclic manner. 



D. Model parameters 

The following quantities are needed in order to solve 
the model equations: the atmospheric density Pfluidj 
gravity grain diameter d and density Pgraim whose val- 
ues are found in the literature and are discussed in the 
next Section; the impact threshold velocity for saltation, 
u*t, and the drag coefficient Cd; the effective restitution 



and 



^sand . 
^0 



7 (eq- 



coefficient a and the heights ^m, 
([T7])) and the wind shear velocity u*. 

The impact threshold velocity is about 80% of the 
threshold for aeolian entrainment, M*ft l7[, which in turn 
is calculated as in Iversen and White Q. This leads to 
the following equation for M*t: 



u*t = 0.8 A. 



Pfluid)fl''^ 



Pfiuid 



(19) 



where A is called the Shields parameter, which depends 
on the shape and sorting of the grains and on the angle of 
internal friction [2l| . The Shields parameter is calculated 
as in Iversen and White ^ : 



A = 0.129 



(1 + 6.0 X 10-Vpg,aingrf' 



0.5 



(l.928ReSr' 
for 0.03 < Re*ft < 10 and 

A 



(20) 



vO.5 



= 0.129(1 + 6.0 X 10-7/9grai„.9d^-^) 
{1 - 0.0858 exp [-0.0617(Re,ft - 10)]} 



(21) 



for Rcttt > 10, where Re*ft is the friction Reynolds num- 
ber Re*ft = Uf:ftd/v, and the constant 6.0 x 10^^ has units 
of kg-m'^'^-s"^, while all other numbers are dimensionless. 
The kinematic viscosity v is defined as ?//pfiuid, where rj 
is the dynamic viscosity. We notice that in contrast to 
V, rj depends only on the atmospheric temperature and 
composition. 

The drag coefficient C'd is a function of the Reynolds 
number Re. Jimenez and Madsen [2^ calculated the drag 
coefficient Cd of a particle falling with settling velocity 
V{ from the balance between the gravitational force and 
the drag resistance of the fluid. To adapt the formula of 



Jimenez and Madsen [23] — which is valid for Re within 
the range 0.2 < Re < 127 — to grain saltation, we con- 
sider the balance between the fluid drag on the grains in 
the saltation layer and the bed friction that compensates 
the grain-born shear stress at the surface. This leads to 
the following equation [13, : 



where 




(/Ogiain - Pfiuid) gd 



Pfluid 



(22) 



(23) 



is called the fluid-sediment parameter and A^ and are 
constants that contain information about the sediment 
shape factor and roundness. In this manner, the drag co- 
efficient for saltating particles is the same as in Jimenez 
and Madsen ^23*] but with the quantity S (eq. ([23]) ) cor- 
rected b y a factor 1/V2a. Furthermore, Jimenez and 
Madsen [23] suggested to use A^ = 0.95 and = 5.12 
for typical applications when particle's shape and round- 
ness are not known. 

The parameters a, and zi are computed using the 
equations obtained by Duran and Herrmann f23] . These 
equations have been obtained from comparison with wind 
tunnel data of Rasmussen et al. [2^ and Iversen and Ras- 
mussen [1^ , and allow to calculate the model parameters 
for saltation in different physical environments. Here we 
just display the equations and refer to the original refer- 
ence [2J| for detailed explanation. 

The equations for the model parameters are scaling 
relations that incorporate the timescale 



and the lengthscale 



I'^Pfiuid 



^^5(/0grain - Pfluid) 



1/3 



(24) 



(25) 



The reference height zi, at which the effective wind ve- 
locity UcS is calculated, is given by the equation 

zi = 35^^. (26) 

The height zi is between the mean saltation height 

= 14u*tii/, (27) 

and the surface roughness 

d/20. (28) 



_sand 



The last equation gives intermediate values between d/30 
[3 and d/8 27]. Finally, the effective restitution coeffi- 
cient a (eq. ([8])) is simply calculated with the formula 

M 



a = 0.17d/C. 



(29) 
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Summary — With the equations presented in this Sec- 
tion, the model parameters can be calculated from the 
following quantities: pfiuid, 9, d, Pgrain, and from the vis- 
cosity 77. Thus, these quantities are, together with 7 (eq. 
(|17p ) and with the shear velocity u■^,, the only parame- 
ters of the model. The threshold velocity u*t is obtained 
with eq. P^. the drag coefficient Cd is given by eq. 



while eqs. 
and a. 



(123) 

are used to obtain zi, 



^sand 



III. SALTATION TRANSPORT ON MARS 

Since the quantities governing saltation are functions 
of the atmospheric conditions, we expect saltation to be 
different depending on the location on Mars. The reason 
is that the average atmospheric pressure and tempera- 
ture may vary within an extremely wide range compared 
to the terrestrial case. In this Section, we estimate the 
average trajectories of saltating grains and the sand flux 
under different atmospheric conditions on Mars. The re- 
sults presented in this Section are then used in the next 
Sections, to calculate formation of dunes on Mars. 



A. Martian atmosphere 

The Mars Global Surveyor Radio Science (MGSRS) 
Team has provided valuable atmospheric data of Mars 
[2^. In particular, the temperature T and the pressure 
P near the surface have been systematically measured in 
many locations of Mars. We use the ideal gas equation to 
calculate the local atmospheric density, pfluid, from the 
MGSRS pressure and temperature data. An atmosphere 
of 100% CO2 is considered. Furthermore, the dynamic 
viscosity 77 of the Martian atmosphere is a function of the 
temperature T, and is calculated using the Sutherland's 
formula (H: 



T + C 



[T/Tq 



,3/2 



(30) 



where for CO2 we have 770 = 1.48-10"^ kg/m-s, C = 240, 
To = 293.15 K j29j]. Finally, the kinematic viscosity 1/ is 
calculated with the equation i/ = r///)fiuid- 



B. Particle size of Martian dunes 

Edgett and Christensen [s^l have used thermal inertia 
data to obtain the grain diameter of dunes in intra-crater 
fields of dark dunes on Mars. They found that the av- 
erage grain diameter of Martian dunes is d — 500 ± 100 
/im, which is coarser than the mean diameter of terres- 
trial dune grains, 250 /im jT9| . The value of mean grain 
diameter d = 500 /im as measured by Edgett and Chris- 
tensen [33 for dunes on Mars is used in the calculations 
of the present work, while d = 250 /xm is considered 



for terrestrial dunes. Furthermore, we take the density 
Pgrain — 3200 and 2650 kg/m"^ for Martian [3l| and ter- 
restrial grains, respectively. 

Why is the sand of Martian dunes coarser than the 
sand of Earth's dunes? There is a critical value of the 
diameter d, below which the particle remains suspended 
in the atmosphere. The critical diameter depends on the 
vertical fluctuating component of the wind velocity u' 
[S^ l- If the standard deviation of u', which scales with 
the wind friction speed, is larger than the settling velocity 
of the grain, V{, then the particle will remain suspended. 
The falling velocity Vf is obtained from the equilibrium 
between the gravitational force and the fluid drag [23j . 
Particles for which the ratio Wf/u*ft is smaller than 1.0 
enter suspension [3^ . In this manner, a critical grain 
diameter of 210 /im is obtained for Mars, while on Earth 
the critical value is about 52 /im (30j . 

The critical diameter obtained in this manner for Mars 
appears inconsistent with the observation of Martian ae- 
olian ripples composed of grains of diameter around 100 
/im [33| . In order to calculate the transition suspen- 
sion/saltation, we use the "falling" velocity Vf obtained 
from the equilibrium between the fluid drag and the bed 
friction Tg. This "falling" velocity is given by the equa- 
tion 




(Pg 



Pfluid)gc? 



Pfluid 



(31) 



in which the value of Cd (eq. ([22])) differs from the drag 
coefficient of a vertically falling grain by a factor of the 
order of l/2a, where a is given by eq. (P^ . 

Figure[2]shows the ratio Vf / u*ft as function of the grain 
diameter d, calculated using parameters for Earth and for 
Mars, where we used the nominal pressure P = 6.0 mbar 
and temperature T = 200 K for Mars. In this figure, 
the threshold shear velocity for saltation, u*ft = 1.25 zi*t, 
is calculated using eq. p^ . while Vf is calculated with 
eq. (pij) . The dashed line indicates the transition value 
V{/u^ft = 1.0. 

As we can see from fig. [5J particles with diameter 
smaller than 45 /im enter suspension on Earth, while the 
critical value of d on Mars is around 110 fim. This value 
is larger than the terrestrial one, but is smaller than the 
one obtained in previous calculations [30| . Furthermore, 
we see that the ratio between the measured average grain 
size of dunes (500 and 250 /im on Mars and on Earth, 
respectively) and the critical diameter obtained in fig. [2] 
is around 5.0 on both planets. In fact, pure saltation 
is expected to occur only if the falling velocity is larger 
than 2.5?i*tj which explains why the sand of dunes is 
effectively much larger than the critical diameter [30| . 



C. Saltation trajectories and sand flux 

The model parameters that govern the grain trajec- 
tories are the average saltation height, Zm (eq. ([77)) '): 
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FIG. 2: Ratio between the falling velocity iif and the thresh- 
old friction speed for saltation, ii.ft, calculated for Mars and 
for Earth as function of the grain diameter d. Dashed line in- 
dicates the saltation/suspension transition at Vi/UfH = 1.0, 
and intercepts the Martian (terrestrial) continuous line at 
d = llO/iin (d = 45^im). 



the reference height zi (eq. (|26p) at which the effective 
wind velocity UcS is calculated; the drag coefficient, Cd 
(eq. (HI])); and the effective restitution coefficient, a (eq. 

From d = 500 /im, wc obtain the surface rough- 
ness zg*^"^ = 25 (eq. (|28l) ). The saturated sand flux 
Qs (eq. (O) is function of the wind shear velocity and 
further depends on the saturation velocity of the saltat- 
ing grains, Vg = \vs\, which is calculated in eq. p2p . The 
quantities controlling saltation on Mars are calculated in 
table H 

In table [H the model parameters have been calculated 
using different values of pressure P and temperature T 
valid for Mars. We see that the minimal friction speed for 
saltation, M,t, on Mars may vary by a factor of 2. We note 
that ranges of P and T even wider than the ones studied 
here may occur on Mars [2^. Moreover, we calculate Vs 
and Qs in table [J using a constant value of w*/w*t = 1-5, 
since this is a representative value for saltation on Earth 
[3^ . The corresponding values calculated for Earth are 
shown in the last row of table HI where we used d — 250 
/im, density pfluid = 1.225 kg/m^ and viscosity 77 = 1.8 
kg/m-s, while g = 9.81 m/s^. 

We see in table |T] that the values of sand flux on Mars 
are typically 10 times larger than on Earth. This is in 
agreement with the findings from wind tunnel simula- 
tions of saltation in a Martian environment by White Q . 
We also see that Martian particles travel with higher av- 
erage velocities, while the mean saltation height Zm on 
Mars is larger than on Earth, and may be over 1.0m de- 
pending on the atmospheric conditions. 



While the ratio u*/u*t in table |T] is constant, in the 
upper inset of fig. [3] we calculate Qs for a constant wind 
velocity = 3.0 m/s using values of P and T within 
the range studied in table |T1 We see that the same wind 
friction speed transports more sand where P is higher 
and T is lower, which means lower pfiuid and u*t. 

In the main plot and the lower inset of fig. [3] we show 
how the flux and the particle velocity at a given location 
on Mars depend on w,/u*t. In the main plot we fix T = 
200 K and calculate the saturated flux Qs for different 
values of atmospheric pressure P as function of M*/u,t — 1- 
In the lower inset, the same calculations are made for 
the grain velocity Vs- We see that the grain velocity in 
equilibrium is determined by the atmospheric conditions 
and has only a weak dependence on the friction speed . 
The equilibrium velocity of the grains, Vs, in fact scales 
with u^f Equation (jl2p can be analytically solved in the 
simple case of the two-dimensional flow over a sand bed, 
where the gravitational term can be disregarded, which 
gives [131 

Vs = Ucff - Vf/V2a. (32) 

Because Woff (eq. pT|) ) and Vf (eq. ([?T|) ) both scale with 
M*t, Vs also does. 

In fig. [3] we see that for a given value of M*/M*t, both 
the flux and the grain velocity are larger for lower atmo- 
spheric pressure P. This is because the shear velocity M*t 
required for sand transport is higher for lower P, while 
Vs scales with and qs scales with w^^ (eq. (HH)). 

Table ini shows Vg and qs calculated for different values 
of u*/u*t on Mars and on Earth. Because the threshold 
wind friction speed on Mars is 10 times higher than on 
Earth, the average velocity of saltating grains on Mars 
is one order of magnitude higher than the velocity of 
Earth's grains. Again, Vg may have different values de- 
pending on the location on Mars, while qs depends fur- 
ther on u^:. 

In summary, using the atmospheric data provided by 
MGS Radio Science Team, we can calculate the quan- 
tities controlling saltation at a given location on Mars, 
for example, at a given dune field. From the "Weather" 
maps 28] , we obtain the value of P and T characteristic 
of the area at which the dune field is located. Next, the 
density and viscosity are calculated from P and T, using 
the ideal gas law and eq. ([50)1 . while the model param- 
eters are obtained, as exemplified in Table IH using the 
grain diameter d = 500 fim of Martian sand dunes. 

The wind velocity in the dune fields on Mars is 
an unknown quantity. It must be determined from the 
calculations of dunes, as we will see in the next Section. 

Indeed, there is still one missing quantity for Mars 
which we need in order to solve the sand transport equa- 
tions: 7, which appears in eq. p6p . 7 is given by the 
product rj, where r (eq. pU)) ) is related to the saltation 
trajectories, and 7 (eq. ([6])) gives the strength of the soil 
erosion. However, r and 7 can not be calculated sepa- 
rately !l3|. It is the quantity 7 (eqs. (fTB|) and HZ])) that 
can be determined from comparison with measurements 
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TABLE I: Main quantities controlling saltation on Mars under several values of pressure P and temperature T, and a constant 
u*/u«t = 1.5. The threshold shear velocity Utt, the mean saltation height Zm, the drag coefficient Cd, and the model variables 
z\ and Q depend on the atmospheric conditions, and have been calculated for a constant grain diameter d = 500 /xm and density 
Pgrain ~ 3200 kg/m^, and with a dynamic viscosity obtained from the temperature (eq. (|30p V The grain velocity Vs and the 
saturated flux have been calculated with eqs. (|32|) and (|14|l . respectively. The corresponding values for terrestrial saltation 
are shown for comparison. On Earth, the value it* = 1.5u,t means a shear velocity of 0.32 m/s. 
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FIG. 3: Main plot: Saturated sand flux as a function of the relative shear velocity u*/utt — 1 for different values of atmospheric 
pressure — and therefore different values of u«t — obtained with a temperature T = 200 K. The lower inset on the right shows 
the corresponding values of the average grain velocity Vs- In the upper inset on the left, we show the saturated flux for u, = 3.0 
m/s calculated for different values of temperature valid on Mars. 



of the transient of flux saturation. The terrestrial value 
7 = 0.2 has been obtained by Sauermann et al. [Tsj from 
comparison of the saturation transient of the flux with 
experimental and numerical data [ssl . IsgI . [37| , which are 
not available for Mars. Therefore, 7 is the parameter of 
the saltation model that remains to be determined for 



saltation transport on Mars. It will be obtained in the 
next Section, from the calculations of Martian barchan 
dunes. 
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TABLE II: Average velocity «s of saltating grains on Earth 
and on Mars as a function of the relative shear velocity u,/ii«t. 
Temperature T — 200 K and pressure P — 6.0 mbar were used 
for Mars. 

IV. DUNE FORMATION ON MARS 

One very common type of dunes on Mars are barchans. 
They have one slip face and two horns, and propagate 
on bedrock under conditions of uni-directional wind (fig. 
|4|. Barchans are the simplest and best known dunes 
0, m, M, M El, 113, El- They are the subject of sci- 
entific and also environmental interest because of their 
high migration speed: on Earth, barchans 1 — 5 m high 
may cover 30 — 100 m in a year. 



windward side 




wind '■■ ^^^^ 

^^^lip face \^ 






■* horns 



FIG. 4: Sketch of a barchan dune showing the windward side, 
horns and slip face. In the inset we see the definitions of dune 
width W and length L. 

On Mars, barchan dunes occur on the floor of craters 
and on the north pole [4?| . Similarly to Earth's barchans, 
they form corridors and large dune fields, as in the 
Arkhangelsky crater (fig. [T^) . It appears surprising that 
intra-crater dunes on Mars look in g eneral similar: they 
have mostly an elongated shape [44| . 

The barchan dunes in the Arkhangelsky crater (41.0°S, 
25.0°W) are amongst the largest barchans on Mars. Fur- 
ther, there are good reasons to begin our study of Mar- 
tian barchans with the Arkhangelsky dunes: they have a 
wide spectrum of dune sizes; have not been significantly 
altered by secondary winds, and do not appear joined at 
their horns forming chains of barchanoids. For example, 
such features are observed in the dunes at Kaiser crater 
and Proctor crater |3l| . 

Let us try to reproduce the shape of the Arkhangelsky 
barchans with the dune model using parameters of the 
present atmosphere of Mars. 

The atmospheric pressure P and temperature T near 
the Arkhangelsky crater are, respectively, 5.5 mbar and 
210 K |28|. These values yield a local Martian atmo- 
spheric density pfluid = 0.014 kg/m"^, and a fluid viscosity 



T] K, 1.06kg/m-s. Using the mean grain diameter d = 500 
/im, grain density Pgrain — 3200 kg/m'^, and gravity 3.71 
m/s^, it follows that the threshold wind friction speed for 
saltation in the Arkhangelsky crater is M,t = 2.12 m/s. In 
this manner, all parameters of the saltation model, but 7 
(eq. (|17p ). are determined using the equations presented 
in the previous Section. 



A. The shape of the Arkhangelsky barchans 

In the calculations of barchans, a constant upwind 
shear velocity > u^t and a small influx qin/Qs (the 
interdune flux) are imposed in x (downwind) direction 
at the inlet, starting with a Gaussian hill having the vol- 
ume similar to the dune we want to reproduce. The hill 
evolves in time until displaying the barchan shape, i.e. 
linear relations between length L, width W and height 

H [i3,lil. 

Measured values of shear velocity of Martian winds are 
mostly between 0.4 and 0.6 m/s [i^, which are much 
lower values than the Martian threshold for saltation 
(« 2.0 m/s). Indeed, many authors estimated that the 
shear velocity of Martian sand-moving winds, which 
occur within gusts of extreme dust storms, may reach 
maximum values between 2.2 and 4.0 m/s 0, @|. Again, 
very unprobably must achieve values of the order of 
4.0 m/s on Mars |5i]. 

On the other hand, we know from experience with ter- 
restrial dune fields that the flux in areas between dunes 
is normally small, between 10 and 40% of the maximum 
flux Qs Moreover, as shown from calculations in pre- 
vious work, the shape of a barchan dune of given size 
depends in an important manner on the interdune flux 
only for values of qin/qs above this range [45|. 

On the basis of the observations obove, we try to repro- 
duce, first, the shape of one Arkhangelsky barchan, which 
has width W ~ 650 m, using q-m/qs = 20%. Further- 
more, we take values of wind friction speed in the maxi- 
mum range between 2.0 and 4.0 m/s, which gives u*/u*t 
approximately between 1.0 and 2.0 in the Arkhangelsky 
crater. Now the quantity 7 remains (eq. (jl7p ). which 
could not be estimated for Mars. As a first guess, we take 
the terrestrial value 7 — 0.2 and investigate whether the 
Arkhangelsky dune can be obtained with in the above 
mentioned range. 

We obtained a surprising result: if we take the same 
7 = 0.2 as on Earth, the Gaussian hill does not evolve 
into a barchan: it simply does not develop a slip face but 
a dome is obtained. However, if we take 7Mars on Mars 
of the order of 10 times the terrestrial value 7Earth = 0.2, 
then barchan dunes with shape similar to the Arkhangel- 
sky barchan can be obtained, as we see in fig. O the elon- 
gated shape characteristic of the Arkhangelsky dunes is a 
result of low values of shear velocity, with w^/it^t smaller 
than 1.5. This means values of up to 3.0 m/s. For a 
constant value of 7 = 2.0, the dune shape deviates from 
the Arkhangelsky barchans for increasing values of u*. 
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FIG. 5; Barchan dunes of width W — 650 m calculated using 
parameters for Mars, with 7 = 10 7Earth, and different values 
of wind shear velocity Ut/utt- 



Thus, the shear velocity in the Arkhangelsky crater must 
be close to the threshold for saltation transport. This ex- 
plains the elongated shape of intra-crater barchan dunes 
on Mars. 



B. Entrainment of saltating grains on Mars 

The wind shear velocity « 3.0 m/s estimated for the 
Arkhangelsky crater is well within predicted maximum 
values of on Mars. But why should the Martian 7 
(eq. (fT7|) ) be ten times larger than on Earth? 

The quantity r = 7/7 (eq. (fTO)l ) on Mars should not 
differ much from the Earth's value. This is because the 
ejection velocity of splashed grains is proportional to the 
velocity of the average impacting grains [46j . which in 
turn scales with the average saltation velocity Vg . In this 
manner, we must understand why the Martian entrain- 
ment rate, 7, differs from the one on Earth. This quantity 
determines the intensity of the grain-bed collisions, the 
modelling of which is beyond the scope of this work [l^ . 

However, Anderson and Haff [i^ showed that the num- 
ber of splashed grains is proportional to the velocity 
Vimp of the impacting grains. Let us rescale Wimp with 
Vcsc — y/gd, which is the velocity necessary to escape 
from the sand bed This velocity has value approx- 
imately 4.5 cm/s, both on Mars and on Earth. Further, 
I'imp scales with the mean grain velocity, Vg. This leads 
to the following scaling relation for the entrainment rate 
of saltating grains: 

7 oc Vs/ \fgd. (33) 

Typical values of the average velocity of saltating grains 
on Mars are shown in fig. [3] and in Table |TT] as function 
of the relative wind friction speed w,/w,t. We see that 
the grain velocity on Mars is one order of magnitude 
larger than on Earth: scales with u,t and has only a 
very weak dependence on which we neglect. In this 
manner, we can write 7 cx u*t/-\/g3. Since we know that 
7 = 0.2 on Earth, where g = 9.81 m/s^, d = 250 /im and 
M*t = 0.218 m/s, we obtain 

7 = 0.0454=- (34) 
^Jgd 

Equation ([M]) gives 7 « 2.24 in the Arkhangelsky crater, 
which is in fact one order of magnitude higher than the 
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FIG. 6: Main plot: characteristic length of flux saturation, 
^s, calculated with eq. (|16p , as function of Ut/u,f The inset 
shows the characteristic time ts ~ is/vs as function of u,/u,t, 
where Vs is the average grain velocity (eq. (|32|l '). 



Earth's value, as obtained in a different way from the 
calculations in fig. [5l 

Summarizing, we found that the entrainment rate of 
grains into saltation on Mars is ten times higher than on 
Earth. This is explained by the Martian larger splash 
events, which are consequence of the higher average ve- 
locity of saltating grains on Mars. 

What is the consequence of a ten times higher entrain- 
ment rate on Mars? Because the saturation length of 
the sand flux depends on the rate at which grains enter 
saltation, the larger splash events on Mars have a crucial 
implication for the formation of sand dunes on the red 
planet. 

While on one hand the lower Martian atmospheric den- 
sity pfluid and gravity g result in longer grain trajectories 
than on Earth , the saturation transient of the flux on 
Mars is shortened by the faster increase in the population 
of saltating grains. This is because the wind strength is 
reduced more rapidly the faster the grains are launched 
into saltation after splash ("feedback effect" [HI). In 
fact, the characteristic length of flux saturation, is (eq. 
(fT6)) ). scales with the average saltation length, £. How- 
ever, is is, furthermore, proportional to I/7. 

The characteristic length of flux saturation, 4, is calcu- 
lated in the main plot of fig. [S]using parameters for Earth 
and for the Arkhangelsky crater on Mars, with u*/u*t in 
the range between 1.0 and 2.0. In this figure, the Mar- 
tian has been calculated using eq. (jl6p with 7 given 
by eq. In the inset of fig. [HI we have calculated 

the characteristic time of flux saturation, ts — (s/vs, for 
different values of u*/u*t. 

It is remarkable that although the Martian and terres- 
trial values of £s differ by a factor of 10, ts on Mars is 
nearly the same as the terrestrial one for a given u*/u*t. 
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FIG. 7: Sand flux calculated over a flat sand bed on Mars 
and on Earth. Main plot: distance of flux saturation, As, 
normalized by the characteristic length 1^, as function of the 
influx Qin/qs at the inlet. Inset: evolution of the normalized 
flux q/q^ with downwind distance x/i^ for gin/gs ~ 0.20. 



This is because the average velocity of saltating grains, 
Us, is one order of magnitude higher on Mars, as shown 
in Section [nil 

As an example, we calculate the sand flux q (eq. 
over a flat sand bed submitted to an unidirectional wind 
of constant strength, using parameters for Earth and 
for the Arkhangelsky crater on Mars. A constant in- 
flux Qin/qs is set at the inlet. The evolution of the nor- 
malized sand flux q/qs with the downwind distance x/£s 
calculated using parameters for Earth (line) and for the 
Arkhangelsky crater (symbols) is shown in the inset of 
fig. [7] with qin/qs — 0.2. In the main plot of fig. [3 we see 
that it takes a distance Ag of about 6 times £s, from the 
edge of the sand bed, for the sand flux to achieve 99% 
of its saturated value qg, using realistic values of qin/qs 
between 0.1 and 0.4. 

If we take, for instance, u^/u^t = 1.50, then we obtain, 
with eq. ((T6)) . £s = 0.71 and 13.62 m on Earth, respec- 
tively in the Arkhangelsky crater. On the basis of fig. 
[3 this leads to As ~ 4.3 m on Earth, while the Martian 
As is approximately 81 m. Furthermore, flux saturation 
is reached within approximately 6ts, i.e. within 3.0 s on 
Earth and 4.6 s on Mars. However, if we had taken the 
terrestrial 7 = 0.2 for Mars, then the values of Martian 
flux transient length and time obtained would be of the 
order of 100, respectively 10 times larger than the terres- 
trial ones. 

The larger value of 7 on Mars shortens the character- 
istic distance of flux saturation (eq. (HH)) by one order 
of magnitude. However, since dunes cannot be smaller 
than the saturation length, this means that the scale of 
dunes that is predicted from the scaling of the flux satu- 
ration distance with the average saltation length, £ (isj . 



is reduced by a factor of 10. In conclusion, the Martian 
larger splash is the missing link to understand the size of 
Martian dunes formed by the thin atmosphere of the red 
planet. 



C. Wind speed and migration velocity of barchans 

In the calculations of fig. [5l wc could estimate the 
shear velocity u^/u^t from comparison with the shape of 
one selected barchan in the field, taking 7Mars = 10 7Earth 
and assuming a given value of the interdune flux qin/qs in 
the Arkhangelsky crater. However, since 7 can be now 
calculated for different atmospheric conditions, we can 
solve the model equations to find the values of u^/u^t 
and qin/qs at a given dune field on Mars. Indeed, both 
field variables can be obtained from comparison with the 
minimal dune [Zsj . as summarized in the next paragraph. 

In order to develop the slip face characteristic of 
barchans, sand hills must reach a minimum size, below 
which they are called domes. As shown in a previous 
work, the minimal dune width Wmin is around 13^s, and 
is approximately independent of the interdune flux, qin/qs 
(isj . In this manner, Wmin yields, through eq. ^6]) . the 
value of u*/u*t at a given dune field. Moreover, once 
u*/u*t is determined, the value of qin/qs can be obtained 
from the shape of the minimal dune: the eccentricity 
p^in/Wmin dccreascs approximately linearly with qin/qs 

In the Arkhangelsky crater on Mars (fig. [1^), the 
minimal dune is indicated by two domes which have 
width VFinin ~ 200 m and length Lmin ~ 400 m. From 
l^min — 200 m, we obtain 4 ~ 15.5 m, which gives 
u*/u*t ~ 1-45 or w 3.07 m/s in the Arkhangelsky 
crater. This is essentially the same result obtained pre- 
viously from comparison with the elongated shape (fig. 
[5]). Next, using this shear velocity, we calculate the eccen- 
tricity Lmin/WVin of the minimal "Arkhangelsky" dune 
as function of qin/qs- We see in fig. [8] that the ratio 
Lmin/Wnun ~ 2.0 is obtained with an average interdune 
flux qin ~ 25% of the saturated flux qg . Again, this value 
is nearly the same interdune flux assumed in the calcu- 
lations of fig. [5] 

Figure [5] shows the results obtained using = 1.45 M*t 
and qin/qs = 0.25. In this figure, we show four Arkhangel- 
sky dunes of different sizes next to dunes calculated with 
the model. Further, the main plot in fig. [9] shows the 
length L as function of width W of the Arkhangelsky 
dunes (circles) and of the dunes obtained in calculations 
(full line). We see that the values of u*/zi*t and qin/qs 
obtained for the Arkhangelsky crater on Mars not only 
reproduce the minimal dune but also describe well the 
dependence of the shape on the dune size. 

In this manner, substituting eq. ((M)) into eq. (|16p . 
we have obtained a closed set of sand transport equa- 
tions that can be solved for different atmospheric condi- 
tions. Furthermore, using the model equations, the value 
of wind friction speed M*/u*t and interdune flux qin/qs in 
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threshold for saltation 0, [i^ . As reported from ob- 
servations of Mars missions, saltation transport on Mars 
occurs during a few seconds in time intervals of several 
years. If, for example, winds on Mars achieve u■^, ~ 3.0 
m/s during 40 s every 5 years Q, then, from fig. [Til we 
see that a Martian barchan of length 200 m would need 
[5 years] -lO'^ ■ (3600 • 24 • 365)/40 « 4, 000 years to move 
1.0 m. This result explains why spacecrafts orbiting Mars 
never revealed any movement of Martian barchan dunes. 
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D. Martian bimodal sand dunes 



FIG. 8: Eccentricity of the minimal dune as function of the 
average interdune fiux, qin/qs, calculated using «*/u«t ~ 1.45, 
P = 5.5 mbar and T = 210 K. We see that the eccentricity 
imin/Wmin ~ 2.0 of the domcs in the Arkhangelsky Crater 
is reproduced with an interdune flux of 25% of the saturated 
flux. 



Barchan field 


Pfluid (kg/m^) 


u,t (m/s) 


Vf, (m/s) 
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0.014 


2.12 


17.8 


2.24 


77.6°N, 103.6°W 


0.022 


1.62 


12.3 


1.71 


Earth 


1.225 


0.22 


1.5 


0.20 



TABLE III: Main quantities controlling saltation on Mars and 
on Earth. 



a given dune field on Mars can be determined from com- 
parison with the shape and the size of the minimal dune, 
on the basis of the results presented in Rcf . 45] . 

Let us study a second Martian barchan field which is 
near the north pole (fig. [TJ)), and where Wniin ~ 80 m 
and imin ~ 130 m. At the location of the field, P = 8.0 
mbar and T — 190 K [2^, and thus u*t ~ 1-62 m/s. 
From the minimal dune width Wmin = 80 m, we obtain 
u*/u*t ~ 1-8 or = 2.92 m/s using eq. (HH). Next, 
in the same way as done for the Arkhangelsky dunes, 
we calculate the eccentricity as function of the interdune 
flux, Qin/qs, and we find that the value imin/W^min ~ 
1.6 is reproduced with Qin/qs = 0.30. In fig. [TUl wc 
see that the behaviour L against W of the barchans in 
this field (stars) is well captured by the model (full line). 
Furthermore, these curves are also shown in the main 
plot of fig. [HI (real and calculated north polar dunes are 
represented by stars, respectively, by the dotted line) for 
comparison with the Arkhangelsky dunes. 

It is interesting that the obtained for the north polar 
field is very similar to that in the Arkhangelsky crater, 
although is lower in the north polar field due to the 
higher pfluid (table |TTT| . 

The values of M*/u*t obtained for Mars are within the 
ran ge o f the ones measured in terrestrial barchan fields 
[13,111, HI] ■ Indeed, we see in fig. [Tl]that, for the same 
value of M*/M*t, Martian barchans would move ten times 
faster than those on Earth. 

However, winds on Mars are only seldom above the 



Because linear and star dunes like the ones found on 
Earth were almost not observed in the first images of 
Mars taken by Mariner 9 and Viking orbiters, it has 
been suggested that non-unimodal wind regimes should 
be very rare on Mars [50] . However, the Mars Global Sur- 
veyor MOC Camera has, more recently, imaged a high 
diversity of dune shapes that had been never observed in 
images of previous missions. 

On bedrock and in areas of low sand availability, there 
appear many exotic and up to now unexplained dune 
forms where barchans should occur if the wind regime 
were uni-directional. Dunes as those in figs. [TJ;— e can- 
not appear in areas of uni-directional winds, for in this 
case barchans should be formed. Indeed, it is possible to 
recognize in the images of figs. [1];— e that the dominant 
winds define a resultant direction. 

We found that the dune shapes in figs. [TJ:— e can be 
obtained with a bimodal wind regime. In our calcula- 
tions, the wind alternates its direction periodically with 
frequency 1/Tw forming an angle 9v, as sketched in fig. 
[12^'. In both directions the strength is the same, namely 
u* = 3.0 m/s, as found from the calculations of barchan 
dunes. In this manner, the value of u*/w*t is particu- 
lar to each field, since M*t depends on the field location 
(table EVl). 

To simulate the change of wind direction, we rotate the 
field by an angle 6'w, keeping the wind direction constant. 
In the calculations, thus, the separation bubble adapts to 
the wind direction after rotation of the field. We use open 
boundaries as in the calculations of barchan dunes. Ini- 
tial condition is a Gaussian hill as before, whose volume 
is taken according to the volume of the dune. 

The angle O^j between the wind directions determines 
which of the different forms in fig. [T^] is obtained. We 
found that a barchan moving in the resulting wind di- 
rection is always obtained if 9^ < 90°. If 9„ « 100°, 
then the dune shape in fig. \T^ ' is achieved. And for 9„ 
of the order of 120° or larger, elongated dune forms as 
those in fig. I12b ^ are obtained, which elongate in time. 
As 6'w — > 180°, a dune form of alternating slip face po- 
sition appears. 

All dune shapes in fig. [12] have been achieved with a 
time Tw in the range of 0.7 to 5.8 days. If the period is 
too large, of the order of a few months, then the dunes 
evolve into barchanoidal forms. 
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FIG. 9: Barchans in the Arkhangelsky crater, 41.0° S, 25.0° W on Mars: Mars Orbiter Camera (MOC) images on the left (image 
courtesy of NASA/JPL/MSSS) and calculated dunes on the right. Main plot: i vs of the Arkhangelsky (circles) and north 
polar barchans at 77.6°N, 103. 6°W (stars). Calculations of Arkhangelsky (north polar) dunes are represented by the continuous 
(dotted) line, obtained with itt/u,t = 1.45 (1.80). The dashed line in the upper inset corresponds to terrestrial dunes obtained 
with u-j^/u-ti = 1.45. In the lower inset, we see L/Lmin vs H/Hmin from calculations of the Arkhangelsky (triangles) and 
terrestrial dunes (dashed line). 
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FIG. 10: Calculations of Martian north polar barchans near 
77.6°N, 103.6°W (fig. [Dd). We see MOC images of dunes 
of different sizes on the left, and on the right we see dunes 
calculated using P = 8.0 mbar, T = 190 K, u. = 2.92 m/s 
and qi-n/qs = 0.29. The plot shows L vs W for the real dunes 
(stars) and for the calculated ones (dotted line), as also dis- 
played in the main plot of fig. [9] for comparison with the 
Arkhangelsky dunes. 
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FIG. 11: Dune velocity Vd as function of dune length L. We 
see that Mars dunes (filled symbols) move typically ten times 
faster than Earth dunes (empty symbols) of same L, obtained 
with similar values of u*/u*t as on Mars. 



The dune shape in fig. [T2k ' has been obtained with 
6»„ 100° and with = 250000 s « 2.9 days, while 
6'w 140° and = 500000 s « 5.8 days has been 
used to calculate the dune in b'. Moreover, we found 
that the structure observed in the dune field of fig. [T2h 
can be obtained by a change in the local wind regime. 
The dune shape in fig. Wlk ' has been obtained in the 



following manner: (i) first, an elongated dune form as 
the one in fig. [T2b ^ is formed with an angle 0w = 120° 
and with = 60000 s w 0.7 day; (ii) next, the angle 
6*^ has been reduced to 80°. Thereafter, the linear dune 
becomes unstable and decays into a string of rounded 
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FIG. 12: MOC images on top (courtesy of 
NASA/JPL/MSSS) and calculations obtained with bi- 
modal wind regimes on bottom. A sketch showing the 
definition of the angle S„ of the wind direction (arrows) is 
shown in a'. The wind changes its direction with frequency 
l/r„. We chose Tw = 2.9 days, 5.8 days and 0.7 day to 
obtain dunes in a', b' and c', respectively. The dune in a' 
has been obtained with 9^ — 100° and the dune in b' with 
= 140°. In c', a linear dune obtained with 9^, — 120° 
decays into barchans after the angle 9^ is reduced to 80°. 
Dune in b' elongates to the right with time, which is not 
observed for the dune in a'. Dune in c' decays further until 
only a string of rounded barchans remains. 



Field 


location 


Pfluid (kg/m^) 


u*t (m/s) 


fig.im 


48.6°S, 25.5°W 


0.017 


1.89 


i,..Lu]. 


49.6°S, 352.9°W 


0.014 


2.06 




76.4°N, 272.9°W 


0.03 


1.35 



TABLE IV: For each dune field in fig. [H the fluid density 
Pfluid and the threshold u,t are calculated from the local pres- 
sure and temperature which are taken from the MGS Radio 
Science data (MGSRS 2006). In spite of the broad range of 
u»t, all dune forms in fig. 1121 have been obtained with one 
single value of u, = 3.0 m/s. 

barchans as seen in fig. [12h . 

It is interesting to notice that our calculations provide 
a different explanation for the formation of the Martian 
dune field in fig. [TJi than that proposed by Bourke [5l| . 
We found that the field in fig. [TJj consists of linear dunes 
which are decaying into barchans, while Bourke [5ll | sug- 
gested an alternative view: the small barchans would 
merge to form the linear dunes. 

The results of fig. [T^] provide evidence for bimodal 
wind regimes on Mars. We find that a variety of Martian 
dune forms which appear in craters and which develop 
on bedrock have been formed by a wind whose direction 
alternates between two main orientations. We conclude 
that if more sand were available in those places, longitu- 
dinal dunes would in fact appear as observed in terrestrial 
sand seas. The study of the shape of linear dunes is the 
subject of a future publication |52l |. 

Again, the wind strength = 3.0 m/s used in the cal- 



culations must be interpreted as the representative value 
of shear velocity that is above the threshold for saltation 
and is responsible for the major changes of the surface 
5]. Because Martian winds are most of the time below 
the threshold for saltation, we expect the timescale Trcai 
of the changes in wind directions to be in reality much 
larger than the values of Tw (a few days) found in the 
calculations. 

We define as /w the fraction of time the wind strength 
is above the threshold for saltation. From the results of 
the calculations of barchans, this means that the Martian 

is around 3.0 m/s a fraction of time. We interpret 
this value of shear velocity as the representative wind 
friction speed associated with the gusts of aeolian activity 
that occur during the strongest dust storms on Mars [J] . 
Further, the real timescale Treai of the changes in wind 
direction is defined through the relation = 2^w/7reai- 

Let us assume that winds above the threshold on Mars 
occur generally during Atgaitation = 40 s at intervals of 
AT = 5 years (2000 days or 1.728 • 10^ s) (|, i.e. 
/w — Aisaitation/AT w 2.31 • 10"''. A characteristic time 
Tw « 1 - 5 days means T„ = 86, 400 - 432, 000 seconds. 
Dividing by Atgaitation — 40 s, this characteristic time 
corresponds to 2,160— 10,800 gusts of saltation trans- 
port. The Martian real time T^eai is 

T 

Treai = -r X 5.0 years « 10, 800 - 54, 000 years. 

^^saltation 

(35) 

Therefore, the real time of oscillation of the wind direc- 
tion on Mars found from our calculations is of the order 
of 10^ years, where it has been assumed that Martian 
saltation occurs as frequently as observed from the Mars 
Missions [l[E[i|. 

V. CONCLUSIONS 

In the present work, we have applied a well established 
dune model, which successfully reproduces the shape of 
terrestrial dunes measured in the field, to study dune 
formation on Mars. In summary, we found that dunes 
observed in the images of Mars could have been formed 
by the action of sand-moving winds that occur occasion- 
ally under the present atmospheric conditions of the red 
planet. Below we give a list of the main conclusions: 

• the quantities controlling Martian saltation, as the 
average grain velocity, mean saltation height and 
saturated fiux may vary in a significant manner de- 
pending on the location on Mars. This is because 
local average values of Martian surface pressure and 
temperature may be very different depending on 
the geographical location; 

• from the shape of barchan dunes on Mars, we found 
that the rate at which Martian grains enter salta- 
tion is 10 times higher than on Earth. The Martian 
higher entrainment rate, which is a result of the 
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larger splash events on Mars, shortens the length 
of flux saturation and reduces the scale of dunes 
that is obtained if only the Martian larger salta- 
tion length is considered [H, ; 

• all dune shapes studied in this work could be repro- 
duced with values of shear velocity that do not ex- 
ceed u» = 3.0 ± 0.1 ni/s, independently of the loca- 
tion on Mars. We interpret this value as the repre- 
sentative friction speed of sand-moving winds that 
occur during the strongest dust storms on Mars; 

• for the same value of relative wind velocity u,/M,t, 
barchans would move ten times faster on Mars than 
on Earth. However, the migration velocity of Mar- 
tian barchans is negligible because saltation trans- 
port in fact occurs only seldom on the present Mars; 

• we found Martian dune shapes that have been 
formed by bimodal wind regimes. The timescale of 
changes in wind direction obtained in calculations 
is of the order of a few days. Taking into account 
that winds transport sand on Mars during some 
tens of seconds in intervals of a few years [y], this 
timescale is in reality of the order of 10, 000—50, 000 
terrestrial years. 

It is interesting to notice that a significant change in 
wind direction (by 90° or more) is expected to occur after 
each extreme of the orbital cycle of Mars, which is deter- 
mined by the combined effect of the precession of its axis 
and the arrival at perihelion [sO, [13, [HE, lEB, [13, [EH] • Be- 
cause of precession of the Martian axis, each pole of Mars 
appears tilted to the sun in a cycle of 51, 000 years. Now 
the latitude which "looks" toward the sun at perihelion 
is 15°S, but this "subsolar latitude at perihelion" (SLP) 
mig rates ±25° about the Equator over a 51kYr time span 
[54| . This orbital cycle is the most important one for the 
climate of Mars, the hemisphere of the SLP having short, 
hot summers, and being the one of the major dust storms 
activity. In 25, 500 years, it is the northern hemisphere 
that will be tilted to the sun. "The large amounts of 
fine dust currently deposited in the northern hemisphere 
in regions such as Tharsis, Arabia, and Elysium will be 
redistributed to the southern hemisphere^^ [59]. This half 
cycle of 25, 500 years is in fact well within the range of 
characteristic time 10, 800 < Trcai < 54, 000 years of bi- 
modal winds found from our calculations of sand dunes 



on Mars. 

In comparison, the timescale Trcai of changes in di- 
rections of bimodal winds in terrestrial dune fields is of 
the order of a few weeks or months. On Earth, linear 
dunes appear due to wind changes that occur seasonally 
[boI . [61| , since the fraction of time /w the wind friction 
speed u* is above the threshold u^t is much larger than 
on Mars [13, [fill- In this manner, we do not expect the 
wind changes due to precession of the Earth's axis to 
play a major role for the shape of terrestrial bimodal 
sand dunes. On the other hand, Martian winds that are 
not associated with the intense dust storms mentioned in 
the last paragraph are too weak to move sand. Indeed, 
such weak winds are responsible for the appearance of 
dust devils that leave ephemeral marks on the surface of 
Martian dunes [3ll |. which appear unmobile. 

The shape of Martian dunes could be only achieved 
with real wind and atmospheric conditions of the present 
Mars because the entrainment rate of grains into salta- 
tion, which we found to be 10 times higher on Mars than 
on Earth, was incorporated in the model equations. In 
fact, it is well known from experiments that the splash 
events on Mars are much larger than on Earth due to 
the higher velocity of Martian grains What we have 
found in the calculations is the implication of the larger 
amount of splashed grains on Mars for the flux saturation 
and formation of dunes. It would be interesting to make 
a full microscopic simulation for the saltation mechanism 
of Mars similar to the one that was recently performed 
by Almeida et al. [63l | to confirm our findings microscop- 
ically. 
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